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ABSTRACT 



Context. 

Aims. Intrinsic absorption is a fundamental physical property to understand the evolution of active galactic nuclei (AGN). Here 
a sample of 1290 AGN, selected in the 2-10 keV band from different flux-limited surveys with very high optical identification 
completeness is studied. 

Methods. The AGN are grouped into two classes, unabsorbed (type-1) and absorbed (type-2), depending on their optical spectroscopic 
classification and X-ray absorption properties, using hardness ratios. Utilizing the optical to X-ray flux ratios, a rough correction for 
the ~ 8% redshift incompleteness still present in the sample is applied. Then the fraction of absorbed sources is determined as a 
function of X-ray luminosity and redshift. 

Results. A strong decrease of the absorbed fraction with X-ray luminosity is found. This can be represented by an almost linear 
decrease from ~ 80% to ~ 20% in the luminosity range log L x =42-46 and is consistent with similar derivations in the optical and 
MIR bands. Several methods are used to study a possible evolution of the absorption fraction. A significant increase of the absorbed 
fraction with redshift is found, which can be described by a power law with a slope ~ (1 + z ~f- 62±0A1 t saturating at a redshift of z~2. A 
simple power law fit ~ (1 + z ) 48±0 08 over the whole redshift iis also marginally consistent with the data. 

Conclusions. The variation of the AGN absorption with luminosity and redshift is described with higher statistical accuracy and 
smaller systematic errors than previous results. The findings have important consequences for the broader context of AGN and galaxy 
co-evolution. Here it is proposed that the cosmic downsizing in the AGN population is due to two different feeding mechanisms: a fast 
process of merger driven accretion at high luminosities and high redshifts versus a slow process of gas accretion from gravitational 
instabilities in galactic disks rebuilding around pre-formed bulges and black holes. 

Key words. AGN absorption - AGN evolution - X-ray background 



1. Introduction 

In recent years it has become obvious that supermassive black 
holes at the centers of galaxies must play an important role in 
the evolution of galaxies. There is a strong correlation between 
the black hole mass and global properties of its host galaxy 
spheroid, like the bulge luminosity (Kormendy & Richstone, 
1995; Magorrian et al., 1998) and the stellar velocity dispersion, 
i.e. the M BH -cr relation (Ferrarese & Merritt 120001 Gebhardt et 
al. 120001 ) . At the same time, evidence is mounting that active 
galactic nuclei (AGN) and galaxies in general undergo very sim- 
ilar evolution patterns. The peaks of AGN activity and star for- 
mation occur in the same redshift range (z = 1.5-2) and there 
is a similar dramatic decline towards low redshift. Moreover, 
the mass density of local dormant supermassive black holes in 
galaxy centers is consistent with the mass density accreted by 
AGN throughout the history of the Universe (Marconi et al. 
2004; Merloni 2004), yielding further evidence for a tight link 
between nuclear black hole activity and the growth of galaxy 
bulges. Many recent theoretical models propose that feedback 
from the growing supermassive black holes plays an important 
role in linking the properties and evolution of central black holes 
and t heir host galaxy (Silk & Rees, 119981 Di Matteo et al., 
2005). The feedback from powerful AGN can quench star for- 
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mation and inhibit the further growth of the most massive galax- 
ies (Scannapieco & Oh, 2004), and this effect is indeed required 
to match the results of hydrodynamical simulations to the ob- 
served bright end of the galaxy luminosity functions (Springel et 
al., [20061 Hopkins et al.,[20Q6). 

A strong dependence of the AGN space density evolution 
on X-ray luminosity has been found, the so called "luminosity- 
dependent density evolution" (LDDE) with a clear increase of 
the peak space density redshift with increasing X-ray luminos- 
ity, both in the soft X-ray ( 0.5-2 keV) and the har d X-ray (2-10 
keV) bands (Miyaji et al. , 120001 La Franca et al„ 120021 Cowie 
et al., [20031 Ueda et al., [20031 Fiore et al., [20031 Hasinger et 
al., |2)05). This "AGN cosmic downsizing" evolution, which re- 
cently has also been confirmed in the radio and optical bands 
(Cirasuolo et al., |2002| Bongiorno et al., |2007l l is very similar 
to the "downsizing" of the faint galaxy population which shows 
a smooth decline of their maximum luminosity with decreas- 
ing redshift for z<l (Cowie et al., 119961 ). The observed "cos- 
mic downsizing" indicates that active star formation and black 
hole growth shift to lower mass galaxies throughout the evolu- 
tion of the Universe, which is somewhat counter-intuitive to the 
standard Cold Dark Matter scenario, where large scale structure 
evolves hierarchically from small to larger entities. The AGN 
results indicate a dramatically different evolutionary picture for 
low-luminosity AGN compared to the high-luminosity QSOs. 
While the rare, high-luminosity objects can form and feed very 
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efficiently, possibly by multiple mergers rather early in the uni- 
verse (see e.g. Li et al., 2006), with their space density declin- 
ing more than two orders of magnitude at redshifts below z=2, 
the bulk of the AGN has to wait much longer to grow or to be 
activated, with a decline of space density by less than a factor 
of 10 below a redshift of one. The late evolution of the low- 
luminosity Seyfert population is very similar to that which is 
required to fit the mid-infrared source counts and background 
(Franceschini et al. 2002]* and also the bulk of the star forma- 
tion in the Universe (Madau 1996 ), while the rapid evolution of 
powerful QSOs traces more closely the history of formation of 
massive spheroids (Franceschini et al. 1999 ). This could indicate 
two different modes of accretion and gas supply for the black 
hole growth with substantially different accretion efficiency (see 
e.g. the models of Cavaliere & Vittorini 2000; Di Matteo et al., 
120011 Merloni |2004l Menci et al., |2004j l. 

A large fraction of the accretion in the Universe is obscured 
by intervening gas and dust clouds. Indeed, the spectrum of the 
X-ray background can very well be described by a combination 
of absorbed and unabsorbed AGN evolving through cosmic time 
(e.g. Comastri et al., [T995l Gilli et al.. 120071 hereafter GCH07). 
However, one of the major uncertainties in the study of AGN 
evolution and black hole growth is the unknown distribution of 
absorbing column densities and its dependence on AGN lumi- 
nosity and on cosmic time. Indeed, studies of local Seyfert 2 
galaxies have shown that a large fraction of these (~40%) is 
Compton thick (Risaliti et al. , 119991 ) and thus practically absent 
in X-ray surveys. Therefore, the luminosity-dependent AGN 
evolution picture could be significantly biased by systematic se- 
lection effects (see e.g. Treister et al., 2004). Ueda et al. (2003 ) 
have for the first time found a significant decrease of the fraction 
of obscured AGN with increasing X-ray luminosity, and similar 
results have been obtained almost simultaneously from indepen- 
dent X-ray selected AGN samples by Steffen et al. (20033 an d 
Hasinger (2004). Recently, these trends have been confirmed in 
optically selected AGN observed in the SDSS (Simpson 120051 
and with Spitzer in the MIR (Maiolino et al., 2007 ; Treister et al. 
2008 ), indicating a break-down of the standard AGN unification 
model (Antonucci 1993 ). The possible evolution of the fraction 
of obscured sources with redshift is still a matter of debate in the 
literature. A tentative redshift dependence of the obscured frac- 
tion was reported by La Franca et al. (2005 ). Other authors (e.g. 
Ueda et al. 120031 GCH07) did not find a sig nificant redshift de- 
pendence. Recently, Treister and Urry (2006 ) claimed a shallow 
increase of the absorbed fraction with redshift. The different re- 
sults could indicate both systematic and statistical effects in the 
analysis. 

In this paper the best available AGN samples selected in 
the 2-10 keV X-ray band with the highest redshift complete- 
ness possible have been compiled, resulting in 1406 objects, of 
which more than 92% have redshift information, mainly from 
optical spectroscopy. This allows for the so far best combination 
of object statistics and completeness in any hard X-ray selected 
AGN sample. In section 2 the hard X-ray selected AGN sam- 
ple is described in detail. Section 3 discusses the classification 
criteria used to discriminate between type-1 (unabsorbed) and 
type-2 (absorbed) AGN. Section 4 presents the number counts 
of the different source populations and section 5 elaborates on a 
method to correct for the small, but significant redshift incom- 
pleteness in the sample. Section 6 is the main body of the paper 
and discusses the luminosity- and redshift dependence of the ab- 
sorbed AGN fraction. Finally, the results are discussed in section 
7 and summarized in section 8. Throughout this work a cosmol- 
ogy with n,„=0.3, Qa=0.7 and Ho=70 km s _I Mpc~' is used. 



2. The hard X-ray selected AGN sample 

For the study of the AGN X-ray absorption and its cosmologi- 
cal evolution, well-defined flux-limited samples of active galac- 
tic nuclei with a high redshift completeness have been selected 
in the 2-10 keV band, with flux limits and survey solid angles 
ranging over five and six orders of magnitude, respectively. 

A total of 1290 X-ray selected AGN were compiled from 
nine independent samples containing a total of 1406 X-ray 
sources selected in the 2-10 keV band. 1106 of these have re- 
liable optical spectroscopic redshift identifications. For the opti- 
cally fainter sources, in particular in those samples with fainter 
X-ray flux limits, it becomes increasingly difficult to obtain re- 
liable spectroscopic redshifts and classifications. However, for 
survey fields with substantial multi-band photometric coverage, 
e.g. the GOODS North and South areas in the CDF-N and CDF- 
S, respectively, it has been shown that photometric redshifts can 
be used reliably (Zheng et al. 120041 Grazian et al. 120061 ) . A total 
of 188 photometric redshifts are used in the samples in Table Q] 
The number of unidentified sources, which do not have a reliable 
redshift determination, either through spectroscopy or through 
photometry, but nevertheless have optical or NIR counterparts, 
is only 112, yielding an identification fraction of 92%. In order 
to control systematic redshift incompleteness effects, crude red- 
shifts have been estimated for the unidentified sources, using the 
X-ray to optical flux ratio, following Fiore et al. (2003 ). The sur- 
veys utilized in this work are summarized in TableQ] Please note 
that the number of type- 1 and type-2 AGN in Table Q] includes 
the photometric and crude redshifts and that other source classes 
like galaxies, stars, clusters etc. are omitted here. 

For several surveys a limit in X-ray flux or off-axis angle 
had to be chosen a posteriori, based on optical completeness 
criteria, thus maximizing the number of sources with redshifts, 
while simultaneously minimizing the number of unidentified ob- 
jects. For the purpose of this paper, all X-ray fluxes in the differ- 
ent samples have been converted to observed 2-10 keV fluxes 
assuming the same spectral index for all sources within each 
sample (see below). While it is in principle possible to attempt 
a correction to intrinsic, rest-frame 2-10 keV fluxes for each 
source individually, using spectral indices and hydrogen column 
densities determined from direct spectral fits or hardness ratios 
for each source, this has not been done here for several reasons. 
This procedure introduces additional statistical errors in the flux 
determination due to the uncertainty in the spectral parameters. 
Most of the sources in the flux-limited samples used here con- 
tain very few X-ray photons so that only a crude hardness ra- 
tio is available that does not allow to determine absorption and 
spectral slope independently, let alone more complicated spec- 
tral models. While these uncertainties have a relatively small 
effect on the flux observed in the source detection band, ex- 
trapolations the source rest-frame would amplify these errors. 
For the comparison with theoretical models it is easier to fold 
these observational effects into the predictions instead of trying 
to attempt to correct the observations. Therefore throughout this 
paper all fluxes and luminosities refer to observed quantities in 
the observed 2-10 keV band. Below we summarize our sample 
selection and completeness for each survey. Figure Q] gives the 
solid angle versus flux curve for the individual surveys and the 
total sample used. 

2.1. The HEAO-1 /Grossan sample 

Shinozaki et al. (2006) have selected a flux-limited sample of 
49 AGN from the HEAO-1 all-sky-survey, which has also been 
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2.2. The ASCA Medium Sensitivity and Large Sky Surveys 
(AMSS/ALSS) 

For the ASCA Large Sky Survey and Medium Sensitivity Survey 
the same subsamples as defined by Ueda et al. (2003) have been 
used in this analysis. The ASCA Large Sky Survey (LSS) is a 
contiguous 7 deg 2 strip in the north Galactic pole region sur- 
veyed with the ASCA GIS and SIS instruments (Ueda et al., 
1999). Optical identifications for 34 sources above a 2-10 keV 
flux limit of 1.2 • 10~ 13 erg cirT 2 s _1 are given by Akiyama et 
al. d2000b and updates in Ueda et al. d2003l . The ASCA Medium 
Sensitivity (GIS) survey is based on the two X-ray catalogues 
published by Ueda et al. ( 12001 1 12005I I. Spectroscopic optical 
identifications are from Akiyama et al. ( 120031 1 for the northern 
subsample and from Ueda (2007, priv. comm.) for the southern 
subsample. In total the ASCA Large Sky and Medium Sensitivity 
Surveys contain 139 objects, of which 125 are AGN and only 
one source remains unidentified. 



Fig. 1. The solid angle as a function of flux limit covered in the 
nine independent surveys utilized in this analysis. 



used in the Ueda et al. (2003) analysis. For the brighter part 
of the sample (28 sources with 2-10 keV fluxes larger than 



2.7 • 10" 11 erg cur 2 s~') they used the classical Piccinotti et al. 
( 1982) sample based on the HEAO-1 A2 instrument on a solid 
angle of 2.7 • 10 4 deg 2 . At fainter fluxes (above 10" erg cm -2 
s~') they augmented this by sources from the catalogue derived 
by Grossan d 1 9921 ) through a crosscorrelation of the HEAO-1 Al 
(Large-Area Counter) and A3 (Modulation Collimator) instru- 
ments (the HEAO-1 A3 MC LASS Catalog). 21 sources from the 
Grossan sample were selected in a sky area of 5105 deg 2 with 
particularly high identification completeness. Shinozaki et al. 
have originally excluded the source H1443+421 (1H1448+415) 
from their low-redshift sample. However, later it turned out that 
the most likely identification for this source is the ROS AT source 
1RXS J144645.8+403510, listed as a QSO at z=0.267 in the 
Simbad data base (Miyaji, 2007, priv. comm.). Including this ob- 
ject, the total number of HEAO-1 AGN used in this analysis is 
50. The fluxes of the Grossan objects have been corrected to first 
order for the bias related to confusion noise by Shinozaki et al.; 
while these authors give fluxes for their sample based on indi- 
vidual spectral slopes, for the analysis presented here, all fluxes 
have been converted assuming a power-law spectrum with a pho- 
ton index of 1.65. 



2.3. The XMM-Newton Hard Bright Serendipitous Sample 
(HBSS) 

The HBSS is a subset of the larger XMM- Newton Bright 
Serendipitous Survey (Delia Ceca et al. 2004) carried out by 
the XMM-Newton Survey Science Center (Watson et al. 12001b 
consortium. It is a complete flux-limited sample of bright X-ray 



-\0' n ergcm.- 1 cts- 
fore 2.15 • 10~ 13 erg cm" 2 s 



sources (F x > 7 ■ 10 erg crrT s ) at high galactic latitude 
(\b\ > 20 deg), selected in the 4.5-7.5 keV energy band. The 
HBSS was designed to have a flat sky coverage of 25.17 deg 2 
at a fixed MOS2 count rate limit of 2 • 10~ 3 cts s~ l in the 4.5-7.5 
keV band. The MOS2 count rates in the HBSS were converted to 
observed fluxes in the 2-10 keV band, assuming a photon power 
law index of 1.9 (using a count rate to flux conversion factor of 
1 . The flux limit of the HBSS is there- 
l . The HBSS sample is now almost 
completely optically identified, leaving only 2 out of 67 objects 
unidentified and thus yielding a spectroscopic completeness of 
~ 97% (Caccianiga et al. 2000, Delia Ceca et al. l2008l . 



2.4. The XMM-Newton Medium-sensitivity Survey (XMS) 

The XMM-Newton Medium-sensitivity Survey (XMS, Barcons 
et al. 120071) is a serendipitous X-ray survey and optical identi- 
fication program of sources with intermediate X-ray fluxes dis- 
covered in 25 XMM-Newton high Galactic latitude fields cov- 
ering a sky area of ~3 deg 2 . For this analysis the XMS-H 
sample, selected in the 2-10 keV band was used. The origi- 
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nal Barcons et al. sample contains 159 sources with 2-10 keV 
fluxes > 3.3 ■ 10~ 14 erg cirT 2 s _I with a spectroscopic identifica- 
tion fraction of 83% (27 unidentified sources). However, the ac- 
tual spectroscopic completeness limit varies from field to field. 
In order to optimise the redshift completeness of the XMS, a 
strategy was chosen, which maximizes the sample of identified 
sources, while minimizing the number of unidentified sources. 
There are 15 fields which either are completely identified down 
to the survey flux limit, or have a small number of unidentified 
sources significantly brighter than the flux limit. For those fields, 
the flux limit was maintained and the unidentified sources enter 
the sample for this analysis. For the remaining 10 fields, uniden- 
tified sources have the lowest flux in the particular subsample. 
In these cases the unidentified sources have been excluded. The 
flux limit for this field was then raised to the geometric mean 
between the flux of the first unidentified and that of the last iden- 
tified source in order to avoid the 'gerrymandering effect' intro- 
duced by defining a post-facto flux limit (see also Hasinger et al. 
2005 ). This effect can be explained in the following way: if one 
has a pre-defined flux limit for a (sub)survey, then the faintest 
source in the sample is expected to lie somewhat above this flux 
limit. If, on the other hand, the flux limit would be defined post 
facto exactly at the flux of the faintest identified source, a small 
bias is introduced. Defining the flux limit between two adjacent 
sources in the original sample minimizes this bias. Table[2]sum- 
marizes the unidentified source statistic for the XMS fields. This 
way a cleaner XMS sample could be defined, comprising 144 
objects and almost half the number of unidentified sources, i.e. 
achieving an identification fraction of 91%. The survey solid an- 
gle has been corrected accordingly (see Fig . [TJ . 

Table 2. Identification completeness in XMS fields used in this 
study 



Number of unidentified Sources 12 3 
Number of Fields 16 6 2 1 



2.5. The Chandra Large Synoptic X-ray Survey (CLASXS) 

CLASXS is a large contiguous survey consisting of nine 
intermediate-depth Chandra pointings in a 3 x 3 grid in the 
Lockman Hole (~ 2.5 deg Northwest of the original LH deep 
survey, see below), in a region with the lowest Galactic neutral 
hydrogen column density (Lockman et al. 1986). A catalogue of 
525 X-ray sources discovered in CLASXS is given in Yang et 
al. (120041 1. Steffen et al. ( 200 4) have presented optical identifica- 
tions and 271 spectroscopic redshifts for the CLASXS, yielding 
an overall spectroscopic completeness of 52%. For the work dis- 
cussed here, the sample of 319 sources detected with a signal- 
to-noise ratio SNR > 2.0 in the hard band (2-8 keV) has been 
utilized. Down to the original survey flux limit of 2.1 • 10~ 15 erg 
cm -2 s _1 , there are 138 objects without redshifts. In order to ob- 
tain an acceptable unidentified fraction, the original catalogue 
had to be cut at a substantially higher flux limit of 1.3 • 10~ 14 erg 
cm -2 s _1 , yielding a sample of 103 sources, of which 20 remain 
unidentified. The redshift completeness is thus 81%. 

2.6. HELLAS2XMM 

HELLAS2XMM (Baldi et al., 120021 is a serendipitous sur- 
vey based on suitable XMM-Newton pointings (complementary 



to the XMS) and building on the experience of the original 
BeppoSAX High Energy Large Area Survey (HELLAS). Fiore 
et al. ( 120031) have presented optical identifications and spectro- 
scopic redshifts for 122 sources selected in the 2-10 keV band in 
five XMM-Newton fields, covering a survey area of 0.9 deg 2 . For 
the discrimination between absorbed and unabsorbed sources 
(see below), hardness ratios from Perola et al. (2004) were used, 
kindly provided by M. Brusa (2007, priv. comm.). Originally 
Fiore et al. published optical spectroscopic redshift identifica- 
tions for 97 of their objects, yielding a redshift completeness of 
~ 80%. Later Maiolino et al. (2006) presented additional red- 
shifts from VLT NIR spectroscopy of optically extremely faint 
objects for two additional sources. Both Fiore et al. and later 
Mignoli et al. (2004) tried to estimate redshifts for the remain- 
ing, optically faint unidentified sources. Fiore et al. employed the 
optical to X-ray flux ratio to obtain crude redshift information. 
Mignoli et al. used further, deep NIR imaging with ISAAC at 
the VLT of optically very faint counterparts in HELLAS2XMM. 
Using both the R-K colours and the morphological information, 
they estimated 'photometric' lower limit redshifts for 9 of the 
originally unidentified sources in the Fiore et al. sample, yielding 
also consistent redshifts for the two objects later identified in the 
Maiolino et al. NIR spectroscopy. In section[5]the X-ray/optical 
flux ratio will be used to estimate crude redshifts for all unidenti- 
fied sources in this study, but for the moment the crude redshifts 
for HELLAS2XMM are not taken into account in the sample 
completeness. However, a 2-10 keV flux limit of 10~ 14 erg cm -2 
s was chosen, thus removing the faintest 4 sources, which are 
all unidentified. 

Recently Cocchia et al. 2007, published photometry and 
spectroscopic redshifts in five additional HELLAS2XMM fields 
providing 59 new redshift identifications for the sample of 110 
new sources. In order to maximize the redshift completeness of 
this sample a procedure similar to the one used for the XMS has 
been applied (see above). However, because in this case some 
fields were observed or identified only outside the central image 
region, an inner and an outer off-axis cut had to be applied in 
addition to a flux cut. As in the case of the XMS, the flux and 
off-axis limits were chosen in the geometric mean between the 
last identified and the first unidentified object. This way flux lim- 
its between 1.18 and 2.56 • 10~ 14 erg cirT 2 s were chosen for 
the different fields and the total additional solid angle collected 
this way was 0.48 deg 2 . The sky coverage of the survey (see 
Cocchia et al. 2007 ) was corrected accordingly. In this way an 
additional 49 sources with only 4 unidentified objects could be 
added to the sample. This leaves in total 171 sources, of which 
27 remain unidentified. The redshift completeness in this sub- 
sample of HELLAS2XMM is thus 84%. 

2.7. The Serendipitous Extragalactic Source Identification 
SEXSI Programme 

SEXSI is a serendipitous survey based on 27 Chandra archival, 
high Galactic latitude pointings followed up spectroscopically 
with the Keck telescopes. Eckart et al. (2006) presented spec- 
troscopic redshift identifications for 438 objects out of a total 
sample of 1034 sources selected in the 2-10 keV band, yield- 
ing an original identification completeness of 42% over a survey 
area of ~ 2 deg 2 . For their sample Eckart et al. have utilized the 
Chandra sources over the whole field-of-view (FOV) and to the 
full sensitivity limit of their observations. In order to optimise 
the completeness of the SEXSI sample to be utilized for this 
study, the original sample was cut both in off-axis angle and in 
flux limit, following a procedure similar to that described above 
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for the XMS. At the outset, the whole sample has been reduced 
to off-axis angles below 8.31', thus concentrating on the inner 
part of the Chandra FOV with higher angular resolution and sen- 
sitivity. Then, a limit was defined separately for each of the 27 
SEXSI pointings both in flux and in off-axis angle in a way to 
limit the fraction of unidentified sources to < 80% in each field. 
A total of 7 out of the original 27 fields had to be rejected, be- 
cause this condition could not be fulfilled. The remaining fields 
had between and 5 unidentified sources left (see Table [3]). In 
order to avoid biases in the effective area versus solid angle 
curve introduced by this 'gerrymandering' procedure, both the 
flux limit and the off-axis limit for each pointing was defined as 
the geometric mean between the corresponding values of the last 
identified and the first unidentified source in each field (sorted in 
flux or in off-axis angle, respectively). In this way, a total num- 
ber of 252 objects were retained for the SEXSI subsample used 
in this study, of which 40 remain unidentified, yielding a spec- 
troscopic completeness of ~ 84% over a solid angle of ~ 1 deg 2 . 

Table 3. Identification completeness in SEXSI fields used in this 
study 



Unidentified Sources 1 2 3 4 5 not used 
Number of Fields 1 8 5 3 2 1 7 



2.8. Deep XMM-Newton survey of the Lockman Hole 
(LH/XMM) 

The Lockmam Hole (LH/XMM) has been observed by XMM- 
Newton 17 times during the PV, AO-1 and AO-2 phases of 
the mission, with total good exposure times in the range 680- 
880 ks in the PN and MOS instruments (see Hasinger et al. 
I200T1 120041 Worsley et al. 120051 and Brunner et al. 120051 for de- 
tails). Spectroscopic optical identifications of the ROSAT sources 
in the LH have been presented by Schmidt et al. (1998} and 
Lehmann et al. ( 120001 ). A catalogue from the XMM-Newton PV 
phase was published by Mainieri et al. ( 120021 ) and some pho- 
tometric redshifts have been discussed in Fadda et al. (2002). 
Here objects were selected from 266 sources detected in the hard 
band of the 637 ksec dataset (Brunner et al. 2008 ) with addi- 
tional spectroscopic identifications obtained with the DEIMOS 
spectrograph on the Keck telescope in spring 2003 and 2004 by 
M. Schmidt and P. Henry (Szokoly al., 120081 ) and also includ- 
ing a number of new photometric redshifts (Mainieri, 2007, priv. 
comm.). In order to maximize the spectroscopic/photometric 
completeness of the sample, objects were selected in two off- 
axis intervals with different 2-10 keV flux limits: 2.5 • 10~ 14 erg 
cm -2 s _1 for off-axis angles in the range 11.0-14.3 arcmin and 
5 • 10~ 15 erg crrT 2 s for off-axis angles smaller than 1 1 arcmin. 
The total number of sources in the LH/XMM survey is 64, with 
4 objects still unidentified. 

2.9. The Chandra Deep Field South (CDF-S) 

For this study the catalogue of Giacconi et al. ( 120021 ), based on 
the 1 Ms observation of the CDF-S (Rosati et al. [2002), was 
used and expanded with two additional objects from the cata- 
logue derived from the same data set by Alexander et al. (2003) 
which were confused in the original Giacconi et al. catalogue, 
but formally fulfill the same detection criteria. One source (XID 
527) was removed, because it was identified with the jet of an ex- 



tended radio source already matched with another X-ray source 
(Mainieri et al., 120081 ). 239 sources detected significantly in the 
2-10 keV band (signal to noise ratio >2.0) and within an off- 
axis angle of 10 arcmin were thus selected. The solid angle area 
curve given in Giacconi et al. (2002) was used, but truncated at 
an off-axis angle of 10 arcmin. 

Spectroscopic identifications in the CDF-S have originally 
been obtained by Szokoly et al. (2004) with the FORS instru- 
ments at the ESO VLT, yielding a spectroscopic completeness 
of ~ 48% in the hard band. Additional spectroscopic redshifts of 
CDF-S X-ray sources have been obtained with the ESO VLT as 
part of the VVDS and GOODS surveys and have been collected 
in the catalogue of Grazian et al. (2006). In the meantime more 
spectroscopic redshifts have been obtained both with VIMOS 
and with FORS2 at the VLT through the GOODS programme 
(PI: C. Cesarsky) and the CDFS follow-up (PI: J. Bergeron). 
Recently also the spectroscopic follow-up observations of the 
Extended Chandra Deep Field South (ECDFS, see Lehmer et 
al., 12000) have yielded additional identifications in the CDF- 
S proper, both from the VLT and from Keck (G. Szokoly, J. 
Silverman, 2007, priv. comm.). For the 239 CDFS sources se- 
lected for this study we have therefore the following break- 
down of spectroscopic redshifts: 111 objects from Szokoly et al. 
(2004 ), 30 redshifts from the public GOODS data releases based 
on VLT FORS and V IMOS spectroscopy (Vanzella et al„ 120051 
Popesso et al., 2008 ), 13 and 4 objects from the ECDFS spec- 
troscopic follow-up using VIMOS at the VLT and DEIMOS at 
the Keck telescopes, respectively (J. Silverman, G. Szokoly et 
al., 2007, priv. comm.), 7 objects from the K20 survey (Mignoli 
et al 120031 ) and 3 redshifts from the VVDS (Le Fevre et al. 
2004). This subsample of the CDF-S has therefore a spectro- 
scopic completeness of 71%. 

The field is also included in the COMBO- 17 intermediate- 
band survey, which gives very reliable photometric redshifts 
for optically brighter sources (Wolf et al., 120041 ). A total of 
9 COMBO- 17 redshifts have been used here. The CDF-S has 
also been surveyed by the HST ACS as part of the GOODS 
(Giavalisco et al. 120051 ), and GEMS (Rix et al. 120051 ) projects. 
Very deep NIR photometry has been obtained with the ISAAC 
camera at the VLT as part of the GOODS project (Mobasher et 
al. 120041 ). Finally, the GOODS area has been covered in the MIR 
spectral range by deep Spitzer observations. The CDF-S there- 
fore offers the highest quality photometric redshifts of faint X- 
ray sources, which have originally been derived by Zheng et al. 
(120041 ) and Mainieri et al. ( 120051 ). More recently, Grazian et al., 
(120061 ) have derived photometric redshifts in the GOODS CDF- 
S area including the Spitzer MIR photometry. Using the better 
infrared photometry, Brusa et al. (2008 ) have obtained new, un- 
ambiguous identifications for two objects, which originally had 
two plausible counterparts in Szokoly et al. (2008): For the X- 
ray source with XID 201 the spectroscopic redshift z=0.679 in 
Szokoly et al. was given for the object 'b' in the Chandra error 
circle, while the Spitzer data indicate that the object 'a' with a 
photometric redshift z p =2.88 (Grazian et al.) is the correct coun- 
terpart. A similar situation exists for the X-ray source with XID 
218, where the spectroscopic redshift z=0.479 for object 'a' has 
been replaced with a photometric redshift z=2.27 for object 'b' 
in the error circle. Several of the original Zheng et al. photomet- 
ric redshifts have been superseded by the better Grazian et al. 
values, however, the overall redshift distribution for the CDF-S 
did not change significantly using the new data. For this analy- 
sis 27 photometric redshifts taken from Grazian et al. were used, 
and 33 redshifts with a photometric quality Q > 0.2 from Zheng 
et al. were utilized. Four objects had too low a quality in their 
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Fig. 2. The hardness ratio of X-ray sources selected from XMM-Newton and Chandra surveys with reliable optical spectroscopic 
classification and hardness ratio errors less equal than 0. 15 as a function of redshift. Left: 320 objects optically classified as broad- 
line AGN (BLAGN). Right: objects without broad lines in their optical spectra. Sources classified as type-1 AGN by their X-ray 
hardness ratios (109 objects) are indicated by filled circles, spectroscopically classified type-2 AGN (200 objects) by open circles 
and objects spectroscopically classified as normal galaxies by green triangles (140 objects). The green solid lines in both diagrams 
are model predictions for a grid of power law AGN spectra with photon index 2, locally absorbed through different column densities, 
as indicated in the plot. 
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Fig. 3. Left: Hard X-ray (2-10 keV) sample of 1290 AGN in the I-Magnitude / 2-10 keV flux plane. The dotted, solid and dashed 
diagonal lines indicate constant X-ray to optical flux ratios of log(fx/f op /=-l,0, and +1, respectively. Filled symbols refer to type-1, 
and open circles to type-2 AGN as defined in Section [3] Right: The same data in the redshift / luminosity plane. Luminosity here 
and throughout the paper refers to the observed luminosity, i.e. uncorrected for intrinsic absorption. 
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photometric redshift solution or no data at all and thus remain 
unidentified, yielding an excellent overall redshift completeness 
of 98%. 

2.10. Chandra Deep Field North (CDF-N) 

A sample of 284 X-ray sources significantly detected in the hard 
band and inside an off-axis angle of 10 arcmin were selected 
from the 2 Ms CDF-N source catalogue by Alexander et al. 
(2003 ). The solid angle area curve given by these authors was 
truncated at 10 arcmin off-axis angle as well. For uniformity with 
the other samples used in this study, the CDF-N fluxes were re- 
computed from the count rates given in Alexander et al., assum- 
ing a photon spectral index of 1.4 for all sources. To convert the 
count rates, which are extracted in the 2-8 keV band, into 2-10 
keV fluxes, a counts to flux conversion factor of 2.88 • 10~ n erg 
cm -2 has been used. 

Optical spectroscopic identifications in the CDF-N have 
been compiled from various catalogues in the literature. The 
largest number of spectroscopic redshifts (1 14) were taken from 
Barger et al. (120011 120031) . and 33 objects from Cohen et al., 
(2000). Another 12 objects were taken from various publica- 
tions summarized in the NASA Extragalactic Database (Bauer 
et al. 120021 Chapman et al. [20031 Cowie et al. [20041 C ristiani 
et al. 2004 Smail et al. [20041 Swinbank et al. [20041 Wirth 
et al. 120041 ) and 6 objects were provided by P. Capak (2007, 
priv. comm.). The redshift of the object VLAJ123642+621331, 
originally placed at z=4.424 by interpreting the single emission 
line in its spectrum as Lya (Waddington et al. 1999, has been 
corrected to z= 1.770, because the photometric redshift of this 
source (Capak 2007, priv. comm.) fits much better, if the line 
is interpreted as [Oil] 3727. With the set of 165 spectroscopic 
redshifts available, this CDF-N subsample has a spectroscopic 
completeness of 58%. Photometric redshifts for the CDF-N have 
been kindly provided for this study by P. Capak (2007, priv. 
comm.). They have been described to some degree in Barger et 
al. d2003 ) and allow this field to be exploited to its full depth. 
Here 111 objects with photometric redshift errors Az<1.5 have 
been used, leaving only 8 objects unidentified. The CDS-N sam- 
ple utilized here therefore has a redshift completeness of 97%, 
very similar to the CDF-S discussed above. 

3. AGN classification criteria 

A crucial prerequisite for the analysis presented here is the dis- 
crimination between absorbed and some controversy in the lit- 
erature. The classical distinction between optically unobscured 
(type-1) and obscured (type-2) AGN is done using a discrim- 
ination through optical spectroscopy. Optical type-1 AGN have 
broad permitted emission lines (>2000 km/s), while optical type- 
2 AGN do not show broad permitted lines, but still have high- 
excitation narrow emission lines. Therefore many works in the 
field are mainly using the presence of broad lines as discrim- 
inator and classify objects as type-1 AGN only if they have 
broad lines (see e.g. Barger et al. 2005; Treister & Urry 2006). 
However, this AGN classification scheme breaks down, when the 
optical spectrum is insufficient to accurately determine the emis- 
sion line widths. At high redshifts and low luminosities there are 
several effects that compromise the optical classification. First, 
in general high redshift objects are faint, so that they require 
very long observing times on large telescopes to obtain spectra 
of sufficient quality allowing to unambiguously discern broad 
emission line components. Secondly, there are redshift ranges, 
where the strong classical broad lines shift out of the observed 
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Fig. 4. Cumulative number counts for different subsamples of 
the data used here. For clarity the data have been normalized to 
a Euclidean slope (N(> S) • S^ 4 ). The upper black band shows 
the data points and lcr errors for the total sample in comparison 
to the analytical description (black dashed line) from Moretti 
et al. (120031 ) . The blue band shows the AGN type-1 number 
counts in comparison to the model prediction (blue dot-short- 
dashed line) by GCH07. The red band is for type-2 AGN, again 
in comparison to the Gilli et al. prediction (red dot-long-dashed 
line). The green band and corresponding data points shows the 
number counts of galaxies (see text), which are consistent with 
a Euclidean behaviour down to the faintest fluxes. The thin solid 
line peaking in the middle of the flux range shows the number 
counts for the ~ 8% unidentified sources in the sample and the 
thin dotted line the counts of those sources with only photomet- 
ric redshifts. 



spectroscopic optical bands. Thirdly, and probably most impor- 
tantly, at high redshifts the spectroscopic slit includes the light 
of the whole host galaxy, which dilutes the spectrum of the AGN 
nucleus. Depending on the ratio between nuclear and host lumi- 
nosity, the host galaxy can easily outshine the AGN nucleus, ren- 
dering the AGN invisible in the optical light. This effect has been 
demonstrated for local bona fide bright Seyfert galaxies, where 
the emission lines are diluted in the host galaxy light when the 
integral light of the galaxy is sent through the spectrograph slit 
(see e.g. Moran et al. 120021 ) and is the main reason, why X-ray 
samples are much more efficient in picking up low-luminosity 
AGN at high redshift. 

The optical spectroscopic AGN classification is mainly 
based on broad permitted lines in the restframe blue and ultravio- 
let part of the spectrum. The amount of obscuration in the optical 
spectrum can be estimated using e.g. the Balmer decrement. At 
Nh > 3 • 10 21 cm 2 the broad UV lines are usually suppressed and 
only narrow lines remain. 

In principle, a classification purely based on X-ray proper- 
ties would be possible. Large hydrogen column densities block 
soft X-rays. Therefore, a suitable AGN classification scheme in 
X-rays could involve the column density A^h that can be deter- 
mined from the X-ray spectra. Usually, faint AGN X-ray spec- 
tra are fit by simple power law models. The effect of increasing 
column density is that the soft part of the spectrum is more and 
more suppressed, i.e. the spectrum becomes harder at high A^h- 
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Table 4. Differential number counts best-fit parameters 



Sample 


N 




7i 




s 6 /io- 4 


72 




x 2 




[deg~ 


2 ] 






[erg cirT 2 s" 1 ] 








Total 


167.9 ± 


15.6 


2.64 ± 


0.03 


0.99 ±0.15 


1.24 ± 


0.06 


1.35 


AGN1 


122.2 ± 


15.0 


2.67 ± 


0.03 


1.39 + 0.22 


0.97 ± 


0.08 


1.36 


AGN2 


35.6 ± 


3.5 


2.62 ± 


0.04 


0.33 ± 0.06 


0.82 ± 


0.15 


1.36 



Many works in the field use an Nh value of 10 to discriminate 
between X-ray absorbed and unabsorbed objects. However, in 
practice it is not possible to determine Nh values for samples 
with typically very faint X-ray sources, due to the small num- 
ber of observed photons. Therefore, one has to resort to hard- 
ness ratios measured from coarse X-ray bands. At high redshifts 
in general it gets more difficult to determine absorption values 
from X-ray spectral information, because the absorption cutoff 
is gradually shifting out of the observed band to lower energies. 
In particular, there are several claims in the literature that at high 
redshifts BLAGN tend to show significant absorption (see e.g. 
Brusa et al.. 120071 Wang et al. . 120071 1. but substantial systematic 
effects have to be taken into account here. First, Nh values can 
only be positive, so that the scatter in the data is producing spuri- 
ous absorption values even for perfectly unobscured objects (see 
e.g. the discussion in the appendix of Tozzi et al.. 120061) . Those 
spurious absorptions get larger with increasing redshift. Also, 
warm, ionized absorbers are present in many AGN, producing 
a deep trough in the X-ray spectrum around the oxygen edge 
(~ 0.5keV), with a continuum recovering at softer energies. If 
the redshift of the object is such that the rest frame soft X-ray 
emission is shifted out of the observed band, these ionized ab- 
sorbers can easily be mistaken as intrinsic cold gas absorption. 
This is in particular true for high-redshift Broad Absorption Line 
(BAL) quasars (see e.g. Hasinger et al., 2004). 

To overcome the difficulties of classification in either the op- 
tical or the X-ray bands, a combined optical/X-ray classification 
scheme has been introduced (Schmidt et al., 1998, Szokoly et 
al. 120041 Zheng et al. |20041 >. It first involves a threshold X-ray 
luminosity, which is set at log(Lx) > 42 erg s _1 to discriminate 
against the most luminous star forming galaxies. All point-like 
X-ray sources with higher luminosities are considered as AGN. 
The second parameter is the X-ray hardness ratio HR, defined as 

HR = (H - S)/(H + S), (1) 

with S and H being the count rates in the soft and hard band, 
respectively. This is illustrated in Figure [2] which shows the 
observed hardness ratios as a function of redshift for objects 
selected from XMM-Newton and Chandra surveys with high- 
quality optical spectroscopy and well-determined hardness ra- 
tios. The differences in the response curves of the different 
XMM-Newton and Chandra surveys are small enough (AHR ~ 
0.1) that they do not significantly distort the observed hard- 
ness ratio distributions. The left part of the figure shows ob- 
jects spectroscopically classified as broad-line AGN (BLAGN), 
while the right part shows AGN without broad emission lines. 
The BLAGN show a very narrow range of hardness ratios, in- 
dependent of redshift. The reason for this is that BLAGN are 
strongly dominated by unabsorbed, powerlaw-type X-ray spec- 
tra, which do not change shape when they are redshifted. Only 
11 out of 320 objects in this sample (3.4%) have HR values 
larger than -0.2, which can therefore conveniently be used as 
a threshold discriminating between X-ray absorbed and unab- 
sorbed objects. BAL QSO are treated as BLAGN in this context. 



They are among the most absorbed BLAGN. On the contrary, 
the non-BLAGN show hardness ratios spread out over a much 
larger range. In particular there is a kind of bimodal distribution, 
about half of the objects have HR values smaller than -0.2, con- 
sistent with unabsorbed spectra, and the rest have larger hard- 
ness ratios. There is also a clear trend of decreasing hardness 
ratios with redshift, which is due to the fact that at higher red- 
shifts the exponential low-energy cutoff in the X-ray spectra of 
absorbed sources is shifting towards lower energies. The latter 
effect leads to a degeneracy between absorbed and unabsorbed 
X-ray sources at high redshifts, which can no longer be discrim- 
inated using their hardness ratios. 

Given this information, both the optical and X-ray selection 
effects and systematic errors in the distinction between absorbed 
and unabsorbed sources can be largely overcome by combining 
optical and X-ray spectral diagnostics. For those sources, which 
are unambiguously classified either as BLAGN or as NLAGN, 
using their optical/NIR spectra, the optical type-l/type-2 clas- 
sification is maintained. For simplicity, Seyfert sub-classes up 
to Syl.5 (i.e. 1.2, 1.4 and 1.5) are sorted into the type-1 class, 
subclasses above the Syl.5 (i.e 1.8, 1.9) are sorted into the type- 
2 class. Sources, for which optical spectroscopy does not allow 
an unambiguous discrimination between type-1 and type-2 are 
classified using their hardness ratios. An AGN with a hardness 
ratio HR < -0.2 is called type-1, and with HR > -0.2 type- 
2. This classification has also the advantage that it can be ap- 
plied for sources with no spectroscopic information, but which 
have reliable photometric redshifts. The other advantage is that 
it is based on a simple measurable quantity, the hardness ratio 
which can in principle be included in models of the X-ray back- 
ground. In addition, there are also obvious systematic effects and 
difficulties. A small fraction of AGN are known not to follow 
the simple optical/X-ray obscuration/absorption correlation, ei- 
ther broad line AGN have absorbed X-ray spectra (e.g. for BAL 
QSOs, see e.g. Hasinger et al., [20021 Chartas et al., [2002b . or 
narrow-line AGN have apparently unabsorbed X-ray spectra (see 
e.g. Panessa & Bassani, 2002). Also, as shown in Figure [2j the 
hardness ratio threshold corresponds to different NH-values as 
a function of redshift. Finally, the dust obscuration affecting the 
optical properties of AGN versus gas absorption affecting the X- 
rays may be different, depending on environment and redshift. 
Nevertheless, these systematic effects are arguably far smaller 
than those of any of the other classifications individually. It is, 
however, important to reiterate that the threshold hardness ratio 
of -0.2 used here is calibrated against the optical discrimination 
between BLAGN and non-BLAGN and therefore corresponds to 
a somewhat lower NH value (~ 3 • 10 21 cirT 2 ) than that typically 
used in the literature. 

Figure [3] shows the AGN sample utilized in this study in 
the optical magnitude-X-ray flux and the redshift-luminosity 
plane, respectively. Type-1 and type-2 objects, classified with 
the above procedure are indicated by filled and open circles. In 
principle it would be possible to calculate emitted rest-frame lu- 
minosities from the known redshift and an absorption column 
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Fig. 5. The optical to X-ray flux ratios of the sample objects as a function of observed X-ray luminosity (left) and X-ray flux (right). 
Objects classified as type-1 AGN are shown by blue open circles, type-2 AGN as filled red circles and galaxies as green triangles. 
Objects with spectroscopic redshift information are plotted with larger symbols and those, for which only photometric redshifts 
exist, with small symbols. The unidentified objects, for which crude redshifts have been assigned, are shown as black asterisks. 



density estimeted from the hardness ratios for each object sep- 
arately. However, the small number of source photons detected 
for most objects in this sample introduces significant statistical 
errors in the results. Due to the positive bound on the N# values, 
systematic biases are introduced into the analysis. Throughout 
this paper, therefore, only observed luminosities in the observed 
2-10 keV band are used. The effects of intrinsic rest-frame lu- 
minosities can always be estimated using a forward modelling 
technique like e.g. in the GCH07 population synthesis model. 



4. Number Counts of different populations 

The combination of a large number of surveys with a wide 
range of sensitivity limits and solid angle coverage presents a 
unique resource. On one hand, the surveys presented here re- 
solve a large fraction of the 2-10 keV X-ray background. On 
the other hand, we have an almost complete optical identifi- 
cation and redshift determination for all components. We are 
thus in the position to study the contribution of different object 
classes to the X-ray background. Using the solid angle versus 
flux limit curve given in Figure Q] the number counts for dif- 
ferent classes of sources were constructed. Figure |4] gives the 
cumulative source counts normalized to a Euclidean behaviour: 
N(>S 2 -io) S^ 4 , where S X i4 is the 2-10 keV flux in units of 10~ 14 
erg cirT 2 s _1 . A Euclidean relation would be a horizontal line in 
this graph. Separate curves are shown for the total sample, the 
type-1 and type-2 AGN, as well as for galaxies. The total and 
the AGN source counts display the well-known deviation from 
a Euclidean behaviour at low fluxes. For each of these popula- 
tions a smoothly broken power was fit to the differential source 
counts: 

n(s) = ( (Ns£ 14 r> + ((N/sf - 7l s% u r l r 1 (2) 



The best-fit parameters for those are given in Table [4] As pre- 
dicted in the GCH07 population synthesis model, type-1 AGN 
dominate at high fluxes, while type-2 AGN start to dominate at 
fluxes below 10~ 14 erg cirT 2 s _1 . There is, however, a significant 
difference between the observed number counts and the model 
prediction, which will be addressed later. 

An earlier attempt to disentangle the contribution of the dif- 
ferent source classes to the number counts has been done by 
Bauer et al.. l2004l (see also Brandt & Hasinger 2005). However, 
there are systematic differences in this analysis: the redshift in- 
completeness and the different classification schemes for ab- 
sorbed and unabsorbed AGN make a direct comparison diffi- 
cult. If we compare the surface densities for different classes 
of sources above a 2-10 keV flux of 2.5 x 10~ 16 erg cm' 2 s' 1 , 
close to the sensitivity limit of the deep Chandra surveys, a 
very good match is found for the total X-ray source population 
(~ 4400 deg' 2 for both analyses). The type-1 AGN surface den- 
sity derived in this analysis (~ 1330 deg' 2 ) is somewhat higher 
than the surface density of X-ray unabsorbed objects in Bauer 
et al. (~ 780 deg' 2 ). On the other hand, the type-2 AGN sur- 
face density here is correspondingly lower (2250 deg' 2 versus 
~ 2900 deg' 2 ), respectively. Figure |4] also shows that galax- 
ies, defined in this analysis as X-ray sources with luminosities 
smaller than 10 42 erg cirT 2 s _1 , appear at the lowest fluxes and 
are consistent with a Euclidean behaviour. Their space density at 
a flux of 2.5xl0~ 16 erg cm' 2 s' 1 is ~ 840 deg' 2 , compared to the 
range of 370 - 640 deg' 2 observed by Bauer et al. The number 
counts of the galaxies defined in this way are about a factor of 
3 higher than those predicted by Ranalli et al. (2003) for purely 
star forming galaxies. This excess could indicate that the objects 
in this sample still have a substantial AGN contribution. Some of 
these X-ray sources could also be associated to diffuse emission 
in groups of galaxies. The figure also shows that the unidenti- 
fied sources are conveniently situated at intermediate fluxes, so 



10 



G. Hasinger: Absorption Properties and Evolution of Active Galactic Nuclei 



that they are unlikely to introduce a bias in the number counts. 
At fainter fluxes a larger proportion of objects are optically too 
faint for reliable spectroscopy, so that photometric redshifts have 
to be employed. 



output both in the optical and X-ray regime, type-2 AGN do in- 
deed show a significant trend of higher X-ray/optical flux ratios 
with increasing redshift. This correlation could be fit by a simple 
proportionality between X-ray luminosity and X-ray to optical 
flux ratio: 
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Fig. 6. Redshift distribution of objects with high-quality spectro- 
scopic redshifts (thick black histogram), with photometric red- 
shifts (thin blue histogram) and with crude redshifts estimated 
using the X-ray to optical flux ratios (dashed red histogram). 



The slope of this relation is the same as in the work of Fiore 
et al. 120031 For all objects without redshift information a crude 
redshift was estimated using this relation, also for the objects 
in the HELLAS2XMM survey, where crude redshifts have al- 
ready been estimated with a similar method by Fiore et al. ( 2003 ) 
and M ignoli et al. ( 120046 . Bauer et al. (12004ft and Barger et al. 
(2005) have shown that the Fiore et al. relation breaks down for 
very faint optical counterparts (R>25). Indeed, a few objects in- 
dicated with asterisks in Figure[5] originating in the CDF-S and 
CDF-N surveys fall off from the main correlation line because 
their optical counterparts are too faint. For these objects, statis- 
tically rather uncertain photometric redhifts have been utilized 
instead of the crude redshifts derived from the Fiore et al. rela- 
tion. 

Figure|6]shows the redshift distribution, separately for spec- 
troscopic, photometric and crude redshifts. It is reassuring that 
the distribution of the crude redshifts is very similar to that of the 
photometric redshifts. Given the fact that the flux distribution of 
the unidentified sources peaks in the middle of the sample range, 
the crude redshifts will likely not bias the results. Although indi- 
vidually they have a very low reliability, in a statistical treatment 
they should give a fair estimate of the systematic effects involved 
in the redshift incompleteness. 

Since the AGN classification scheme discussed above, in the 
absence of optical spectroscopy can still discriminate between 
type-1 and type-2 AGN, it is also applied to the unidentified 
sources with crude redshifts in order to estimate luminosities 
and to study the effect of the incompleteness on the analysis dis- 
cussed below. 



5. Correcting redshift incompleteness 

As noted in Table Q] about 8% of the sample remains without 
redshift identification, compared to about 13% with only photo- 
metric redshifts. While this is a small fraction, the incomplete- 
ness most likely depends on redshift, so that systematic redshift 
and therefore spurious evolution effects may be introduced in 
the analysis. On the other hand, as Figure [4] shows, the uniden- 
tified sources are mainly at intermediate fluxes. The reason is 
that they are typically from serendipitous Chandra and XMM- 
Newton surveys, where the quality of the optical photometry 
available does not allow a photometric redshift determination as 
accurately as in the deepest fields, which have excellent multi- 
band photometry. 

It has been shown by Fiore et al. 120031 that for objects, 
where the optical light is dominated by their host galaxy (i.e. 
non Broad-line AGN), there exists a correlation between X- 
ray luminosity and the optical to X-ray flux ratio f x /f pt- This 
is probably mainly due to the fact that AGN host galaxies have 
a rather limited range in absolute luminosity, while AGN lumi- 
nosities can occupy a much larger range. In order to test this 
conjecture for the AGN sample utilized here, the f x /fop t ratios are 
plotted against luminosity and X-ray flux in Figure[3] As is well 
known, the X-ray sources show a substantial scatter in their X- 
yay to optical properties. However, unlike type-1 AGN, which 
at high luminosity are dominated by their non-thermal energy 



6. The fraction of absorbed AGN 

6. 1 . Luminosity-dependence of type-2 fraction 

A comparison of the open and closed circles (AGN type-2 and 
type-1, respective ly) in the right panel of figure [3] shows that 
there is a much larger fr action of type-2 AGN at lower luminosi- 
ties. Indeed, this confirms what had been found previously by 
different teams. Ueda et al. (2003) studied the h ard X-ray lumi- 
nosity functions of ~230 AGN, based largely on ASCA and deep 
Chandra s urveys, with about 95% completeness. They found 
that the fraction of type-2 AGN decreases with X -ray luminos- 
ity, L2-iokeV- Similar results have been obtained independently 
from diffe rent samples by Steffen et al. ({2003 > and Hasinger 
(120041 and later by La Franca et al. (120051) . 

To demonstrate the trend observed in the sample studied 
here, the type-2 fraction (i.e. the number of type-2 AGN divided 
by the total number of AGN in any particular subsample) was 
first calculated as a function of X-ray luminosity. In order to 
minimize possible redshift effects, the analysis was carried out 
in the redshift range 0.2 < z < 3.2. This redshift range was cho- 
sen for simplicity to be compatible with the redshift shells used 
in the later analysis and to minimize any systematic effects that 
could happen at the lowest and highest redshifts. The results are 
shown in Tableland Figure|7](left). In order to keep track of sys- 
tematic selection effects introduced by the unidentified sources, 
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the calculation was performed separately for the sample ignor- 
ing the unidentified sources (thin blue error bars in the figure and 
quantities with primes in the table) and for that with crude red- 
shifts included (thick black symbols in the figure and quantities 
without primes in the table). Because the individual subsamples 
of type-1 and type-2 AGN sometimes contain small numbers, 
for which Gaussian statistics do not apply, the Gaussian equiv- 
alent error bars of the type-2 fraction were computed following 
Gehrels dl986b : AN = + 0.75. 



Table 5. Type-2 fraction as a function of luminosity, measured 
in the redshift interval z=0.2-3.2 



Lmin 


Lmax 


tl' 


t2' 


t2'/(tl'+t2') 


tl 


t2 


t2/(tl+t2) 


42.0 


42.2 


4 


18 


0.82+0.08 


5 


18 


0.78+0.08 


42.2 


42.5 


8 


40 


0.83±0.05 


9 


40 


0.82+0.06 


42.5 


42.7 


24 


38 


0.61+0.06 


24 


39 


0.63+0.06 


42.7 


43.0 


25 


52 


0.68+0.05 


26 


53 


0.67+0.05 


43.0 


43.2 


38 


40 


0.51+0.06 


41 


42 


0.51+0.06 


43.2 


43.5 


48 


52 


0.52+0.05 


54 


53 


0.50+0.05 


43.5 


43.7 


58 


65 


0.53+0.05 


65 


72 


0.53+0.04 


43.7 


44.0 


71 


54 


0.43+0.04 


71 


64 


0.47+0.04 


44.0 


44.2 


93 


31 


0.25+0.04 


103 


41 


0.28+0.04 


44.2 


44.5 


121 


20 


0.14+0.03 


128 


34 


0.21+0.03 


44.5 


44.7 


88 


14 


0.14+0.03 


92 


21 


0.19+0.04 


44.7 


45.0 


58 


12 


0.17+0.05 


58 


17 


0.23+0.05 


45.0 


45.2 


25 


3 


0.11+0.06 


26 


6 


0.19+0.07 


45.2 


45.5 


16 


1 


0.06+0.06 


16 


1 


0.06+0.06 


45.5 


46.2 


13 


1 


0.07+0.07 


13 


1 


0.07+0.07 



Quantities tl' and t2' with primes refer to the sample ignoring the 
unidentified sources, while quantities tl, t2 without primes refer to the 
sample including crude redshifts for the unidentified sources. 



Figure [7] (left) shows the trend of a decreasing absorp- 
tion fraction with increasing X-ray luminosity with unprece- 
dented accuracy. This figure also demonstrates that the red- 
shift incompleteness effects are moderate and mainly affect ob- 
jects at higher luminosity. A simple linear function has been fit 
to the completeness-corrected data points. It yields a slope of 
-0.226 ± 0.014 and a marginally acceptable reduced^- 2 of 1.96. 
A flattening of the decline is anyhow expected at the highest lu- 
minosities, because the simple linear fit does not preclude nega- 
tive values. A comparison of this new data with other results in 
the literature will be made in SectionfT] 

6.2. Redshift-dependence of AGN type-2 fraction 

In order to address a possible evolution of the obscuration frac- 
tion with redshift, one can first simply determine the observed 
ratio of type-2 versus total AGN, integrated over all luminosi- 
ties, in shells of increasing redshift. This is very similar to what 
was done by Treister & Urry |2006). Figure [7] (right) shows 
the dependence of the observed fraction of absorbed sources on 
redshift. To first order this data indicates a flat behaviour, apart 
from a small increase at redshifts below z=0.8. The diagram also 
shows the type-2 fraction as a function of redshift predicted for 
the current sample from the most recent background synthesis 
model (GCH07), which assumes no evolution in the absorbed 
fraction. This curve shows that the small rise of the observed 
type-2 fraction at the lowest redshifts can be understood as an 
effect of the different flux limits for type-1 and type-2 sources in 
each survey (the latter ones are harder to find because of absorp- 



tion). The comparison with this model and with the Treister & 
Urry data points will be discussed in Section|7] 



Table 6. Type-2 fraction norm in redshift intervals 



z Range 


Norm 


corr 


corr. Norm 


Slope 


0.015-0.2 


0.224 ± 


0.045 


0.83 


0.269 


± 0.054 


-0.25 ± 


0.06 


0.2-0.4 


0.273 ± 


0.041 


0.87 


0.313 


± 0.047 


-0.16 ± 


0.05 


0.4-0.8 


0.390 ± 


0.025 


0.95 


0.409 


± 0.026 


-0.28 ± 


0.03 


0.8-1.2 


0.394 ± 


0.027 


1.00 


0.396 


± 0.027 


-0.29 ± 


0.03 


1.2-1.6 


0.468 ± 


0.030 


0.98 


0.476 


±0.031 


-0.30 ± 


0.04 


1.6-2.4 


0.564 ± 


0.030 


1.00 


0.565 


± 0.030 


-0.36 ± 


0.05 


2.4-3.2 


0.575 ± 


0.060 


1.04 


0.551 


± 0.058 


-0.27 ± 


0.11 


3.2-5.2 


0.615 ± 


0.086 


1.06 


0.578 


±0.081 


-0.11 ± 


0.12 



Figure|7](right) is integrated over all luminosities in each red- 
shift shell and is therefore possibly hiding trends, which might 
be observable in a better resolved parameter space. For a bet- 
ter diagnostic of the redshift evolution of the type-2 fraction, the 
dependence of the type-2 fraction on luminosity was analysed 
separately in different redshift shells. In Figure[8]the same trend 
of a decreasing absorption fraction as a function of X-ray lu- 
minosity, as observed in the total sample (see Figure [7] left) is 
confirmed in each individual redshift shell. The linear fit to the 
total sample in the redshift range 0.2-3.2 from Figure [7] left) is 
shown as a dotted blue line in each panel. The same linear trend 
has been fit to the data points in each individual panel. First both 
the normalization of the curve at log(Lx) = 43.75 as well as the 
slope of the relation have been left as free parameters in the indi- 
vidual fits. The last column of Table [6] shows the slopes derived 
for the individual fits. Since they are all consistent with each 
other, the average slope over all panels has been determined to 
-0.281 ±0.016. s slope is somewhat steeper than that derived for 
the total sample in the redshift range 0.2-3.2 (see above), which 
is a consequence of the increasing absorbed fraction with red- 
shift. 

The data in the individual panels have then been again fit 
with a linear relation, but this time keeping the slope fixed to the 
average value of -0.28 1 . The results of these fits are shown as the 
solid red lines in Figure[8]and also in Table|6] For every redshift 
shell this table shows in column (2) the norm of the linear fit to 
the luminosity-dependence, taken at a value of log(Lx) = 43.75 
in the middle of the range. There is a clear trend of an increasing 
normalization of the absorbed fraction with redshift which can 
be seen in the comparison of the blue dotted and red solid lines 
in Figure [8] and in the corresponding data in the table. However, 
in order to quantitatively estimate the evolution of the absorbed 
fraction with redshift, the systematic selection effects have to be 
corrected for. They are due to the fact that in every flux-limited 
sample it is harder to detect absorbed objects compared to un- 
absorbed ones of the same intrinsic flux. Folding this selection 
effect over the solid angle as a function of flux limit of the cur- 
rent meta-sample, one can in principle obtain a correction curve 
as a function of redshift. For simplicity, here the green curve in 
Figure [7] (right) determined for the GCH07 model folded with 
the effective solid angle sensitivity curve in Figure [TJ and nor- 
malized to a redshift of 2 has been used. The (small) correction 
values are given in column (3) of Table [6] With this correction 
applied, finally the corrected normalization can be calculated for 
every redshift shell, wich is given as column (4) of Table|6] 

This corrected normalization has been plotted as a func- 
tion of redshift in Figure [9] A significant increase of the ab- 
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Fig. 7. Left: Type-2 fraction as a function of luminosity, based on the data in the redshift range z=0.2-3.2. The thin blue data points 
show the sample ignoring the unidentified sources, while the thick black data points with filled circles include the crude redshifts 
derived in Section [5] For clarity the thin blue points are shifted slightly to higher luminosity. A simple linear function fit to the 
black data points is shown as a dashed black line. The data points with red circles are derived using optical line widths of a sample 
of [OIII] selected AGN from the SDSS (Simpson 2005). The solid magenta line is from silicate dust studies of a sample of AGN 
observed with Spitzer in the MIR (Maiolino et al. 120071 1. The dotted green line shows the ratio of Compton-thin absorbed AGN 
to all Compton-thin AGN assumed in the GCH07 population synthesis model. Right: The observed fraction of type-2 AGN as a 
function of redshift. The colour coding of the thin blue and thick black data points is the same as in the left figure. The green solid 
line shows the prediction of the GCH07 model, assuming no redshift evolution in the absorbed fraction and folding over the solid 
angle sensitivity curve from FigureQ] The error bars along the Y-axis in this and later figures all give l-cr uncertainties. 



sorbed fraction norm with redshift is seen from this figure. This 
increase, however, seems to saturate above a redshift of z~2. 
Therefore, a broken power law has been fit to the data: t 2 {z) = 
f 2 (0) • (1 + z) ai for z < Zb , saturating at t 2 (z) = t 2 (z b ) for z > Zb- 
The best-fit parameters are determined to ?2(0) = 0.279 + 0.024, 
a\ = 0.62 + 0. 11, Zb = 2.06 + 0.47 and the reduced^ 2 ofthisfitis 
0.82. The evolution in the absorbed fraction in the redshift range 
up to z~2 is therefore confirmed with a statistical significance 
of 5.7<r. In order to assess the significance of the apparent evo- 
lution saturation at high redshifts and to compare with previous 
works, also a simple power law t 2 (z) = t 2 {0) • (1 + z) ai has been 
fit to the whole redshift range. The best-fit parameters for this 
representation are f 2 (0) = 0.308 + 0.022 and a 2 = 0.48 + 0.08 
(i.e. a significance of 6.1cr). With a reduced^- 2 of 1.30 this fit is 
worse than for the broken power law model but still marginally 
acceptable. Finally, a constant value can be ruled out with very 
high confidence (reduced^ 2 = 6.37). 

A simple power law t 2 (z) - t 2 {0) ■ (1 + z)" 3 has also been fit 
to the data without including the crude redshifts (thin blue data 
points in all diagrams). The best-fit parameters for this data set 
are t 2 (0) = 0.309 + 0.024 and a 3 = 0.38 + 0.09 (a significance 
of 4.4<x) with a reduced ;f 2 of 1.08. Including the crude redshift 
thus pronounces the evolution in the redshift range up to z=2 
and provides evidence for a saturation of this evolution. In order 
to assess the significance of this saturation it is interesting to 
note the actual numbers of objects in the last two redshift bins, 
both including and excluding the crude redshifts. The redshift 



bin 2.4 < z < 3.2 contains 39 (38) type-1 and 29 (26) type- 
2 AGN (numbers in brackets refer to the sample without crude 
redshifts). The redshift bin 3.2 < z < 5.2 contains 23 (19) type-1 
AGN and 13 (7) type-2 AGN. The inclusion of crude redshifts is 
thus affecting mostly the highest redshif bin. 

7. Discussion 

7. 1 . Systematic effects of different AGN classification 
schemes 

In section[3]different classification techniques for X-ray selected 
AGN have been discussed. In most classical treatments, type-1 
AGN are just defined by the presence of broad permitted lines 
observed in their optical spectra (see e.g. Steffen et al., 2003; 
Barger et al.,|2005; Treister & Urry |2006j l. In order to estimate 
the systematic differences with the combined optical/X-ray clas- 
sification employed in this paper, the luminosity- and redshift 
dependence of the type-2 fraction has been recomputed for the 
sample using the same BLAGN classification scheme. In this 
case type-1 AGN are spectroscopically identified objects with 
broad emission lines, while type-2 AGN are all other spectro- 
scopically identified objects. For consistency with other works, 
objects with only photometric or crude redshifts are ignored 
here. Figure [10] shows the results of this analysis in compari- 
son to the combined X-ray/optical classification and to the work 
from Treister & Urry (2006) and Treister et al. (2008), respec- 
tively. The left panel shows the observed dependence of the type- 
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Fig. 8. Variation of the type 2 fraction with luminosity for different redshift shells. The thin blue data points show the sample 
ignoring the unidentified sources, while the thick black data points with filled circles include the crude redshifts derived in Section 
[5] The dotted blue line gives the fit to the average relation in the redshift range 0.2-3.2. 



2 fraction on X-ray luminosity. The data points derived from the 
BLAGN only classification (green) are fully consistent with the 
data from Treister et al. (2008), but at lower X-ray luminosi- 
ties both are significantly higher than those from the combined 
X-ray/optical classification used here. This is mainly due to the 
large number of unabsorbed AGN at low to intermediate lumi- 
nosities, missing in the BLAGN sample presumably due to dilu- 
tion from the host galaxy light (see Figure|2|. The same trend is 
also seen in the right panel of Figure [TUJ which shows a signif- 
icant difference in the redshift dependence of the the absorbed 
fractions determined by the two different methods. 



7.2. Luminosity dependence of absorbed fraction 

All analyses represented in Figure [7] (left), based on X-ray and 
optical samples, agree that there is a strong decline of the ab- 
sorbed fraction with AGN luminosity. The X-ray analysis pre- 
sented here confirms this trend in much more detail and with 
better statistical and systematic errors than previous X-ray anal- 
yses. Similar trends are found in other wavebands. Local Seyfert 
galaxies selected from the SDSS spectroscopy sample using the 
[OIII] lines show a fraction of BLAGN increasing with the lumi- 
nosity of the isotropically emitted [O III] narrow emission line 
(Simpson 2005). In order to convert the [O III] line luminosity 
into a 2-10 keV luminosity for the display in Figure [7] the em- 



pirical luminosity-dependent correlations found by Netzer et al. 
(120061) for type-1 and type-2 Seyfert galaxies were averaged: 

log = 13.7 - 0.358 logL 2 _ 10 . (4) 

1-2-10 

Mid-infrared (MIR) Spitzer low resolution spectra of 25 lu- 
minous QSOS at redshifts between 2 and 3.5, combined with 
low-luminosity type-1 AGN confirm this picture (Maiolino et al. 
120071) . These authors find a strong, non-linear decrease (a factor 
of ten) of the thermal continuum dust emission at 6.7/im with 
increasing 5100 A optical continuum emission and [OIII]5007 
line luminosity. They interprete this as a decreasing covering 
factor of circumnuclear dust as a function of luminosity. They 
clearly detect silicate dust emission features (see also Sturm et 
al., [2006), whose strength correlates with luminosity, accretion 
rate L/Lgda and black hole mass Mbh- 

Independently, Treister et al. (2008 ) recently attempted to 
measure this covering fraction by studying the ratio of the MIR 
to bolometric luminosity in a sample of z~l BLAGN selected 
optically in the Sloan Digital Sky Survey (SDSS), the Great 
Observatories Origins Deep Survey (GOODS) and the Cosmic 
Evolution Survey (COSMOS). They use the Spitzer 24/j flux as 
a proxy of the MIR dust-reprocessed radiation. For their redshift 
range this corresponds to ~ 12/i, which according to Maiolino 
et al. ( 120071) and Sturm et al. (2005 ) is mostly dominated by the 
silicate emission and less due to the reprocessed thermal dust 
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Fig. 10. Comparison of type-2 fractions as a function of luminosity (left) and a function of redshift (right) for different AGN 
classification schemes. The filled black data points with thick error bars correspond to the combined optical/X-ray classification 
introduced in section[3]and are identical to those in Figure[7j The green data points with open squares and thin error bars have been 
derived using a purely optical spectroscopic classification, where only BLAGN are classified as type-1 objects (see Figure|2]i. Red 
data points with dotted error bars are from Treister et al. (2008) (left) and from Treister & Urry 2006 (right), respectively. 



emission. The silicate emission is thought to come mostly from 
the narrow-line region (NLR) and thus does not directly probe 
the covering factor of the molecular dust torus. To estimate the 
bolometric luminosity of their AGN (primarily the blue bump) 
they use the Galex NUV data, sampling the rest frame waveband 
875-1400A, which is possibley affected by absorption from the 
Ly-forest and even small amounts of dust extinction. Given these 
systematic differences, it is surprising that their results almost 
perfectly match the Maiolino et al. curve in Figure [7J (left). 

Of course, all the sample selections represented in Figure [7J 
(left) are subject to systematic errors and selection biases. X- 
ray surveys in the 2-10 keV band almost completely miss the 
population of Compton-thick Seyfert galaxies, which are barely 
detectable only above 10 keV (see e.g. GCH07). However, first 
significant minority samples of Compton-thick candidate ob- 
jects have been detected in deep-wide X-ray surveys (Tozzi et 
al., 120061 Hasinger et al., [20071 Brunner et al., |20U8l , consis- 
tent with the GCH07 model predictions. Including information 
from deep MIR surveys, there is now convincing evidence that 
the population of Compton-thick intermediate-luminosity AGN 
even at higher redshifts is of the same order as that of Compton- 
thin AGN (see e.g. the recent results in the CDFS by Daddi et al. 
120071 and Fipreet al. |2008a| as well as in the COSMOS field by 
Fiore et al. 2008b), again fully consistent with the predictions of 
the GCH07 model. Optical surveys, on the other hand, while in- 
cluding many of the Compton-thick AGN missing in X-ray sam- 
ples, fail on the AGN populations diluted by their host galaxy 
(see Figure [2]), and also those so heavily obscured that very lit- 
tle scattering or ionizing radiation escapes (see e.g. Ueda et al., 
120071 1. Nevertheless, the agreement of better than 20% between 
the optical sample selected from the SDSS using the [OIII] lines 



(Simpson 2005) and the hard X-ray selection from this paper is 
reassuring in this respect. 

In Figure [7J also the relation between absorbed Compton- 
thin AGN to all Compton-thin AGN assumed in the GCH07 X- 
ray background population synthesis model is shown. Because 
Compton-thick objects have been explicitely excluded in this 
curve, it can be directly compared to the observed 2-10 keV 
X-ray sample results. Please note, however that this is the as- 
sumed intrinsic ratio, the real observations are still affected by 
additional selection effects. Nevertheless, it is apparent that this 
model overpredicts the fraction of obscured sources, in partic- 
ular at high luminosities. This is also likely one of the reasons 
for the systematic difference between the predicted and observed 
number counts for type-1 and type-2 AGN in Figure[4] Probably 
the best approximation of the "true" fraction of obscured objects 
is the Maiolino et al. relation. Hopefully, the GCH07 model can 
be updated accordingly in the future. 

Physically, this luminosity-dependent fraction might be in- 
terpreted as a 'cleanout effect' because more luminous AGN can 
dissociate, ionize and finally blow away the dust in their envi- 
ronment. At face value this trend signals a break-down of the 
simple unified AGN model, where the difference between type- 
1 and type-2 AGN is solely in the geometry of the observer's 
line of sight with respect to the absorber (Antonucci 1993). The 
luminosity dependence of the absorbed fraction indicates an in- 
trinsic physical difference in the average absorber properties as 
a function of intrinsic AGN luminosity. The simplest interpreta- 
tion of these findings is that low-luminosity AGN are surrounded 
by an absorbing/obscuring medium which covers a large solid 
angle (~ 80% of the sky as seen from the black hole), however 
with an average column density much lower (Nh ~ 10 22 cm~ 2 ) 
than the putative Compton-thick dusty torus of the unified model 
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Fig. 9. Dependence of the normalization of the type-2 fraction 
at log(L x )=43.75 with redshift (from the fits shown in Figure[8]). 
The red solid curve shows a saturated power law fit to the solid 
black data (see text). The blue dotted lines give the 90% confi- 
dence intervals for this fit. The dashed red line shows a single 
poer law fit to the black solid data over the full redshift range, 
while the dot-dashed line shows the same single power law fir to 
the thin blue data points without including the crude redshifts. 
The dashed green data points are from La Franca et al. d20051 >. 



(N H ~ 10 24 - 25 cm- 2 ). High-luminosity AGN are then able to 
clean out their environment, either by ionizing the surround- 
ing medium, or by blowing it away through an outflowing wind 
(see also Miiller & Hasinger 2007). This behaviour is also a key 
ingredient in the "blast wave AGN feedback" model presented 
very recently by Menci et al. (2008), where AGN accretion lu- 
minosity, feedback and absorption are intimately connected to 
each other. 



7.3. Redshift dependence of absorbed fraction 

The evolution of the luminosity-dependent obscuration fraction 
with redshift is still a matter of debate. A possible redshift depen- 
dence of the obscured fraction was reported by La Franca et al. 
(120051 1. using the HELLAS2XMM sample combined with other 
published catalogs, based on a total of 508 AGNs selected in the 
2-10 keV band. On the other hand, other authors (e.g. Ueda et 
al. 120031 GCH07) did not find a significant redshift dependence. 
Recently, Treister and Urry (2006) performed an analysis, using 
an AGN meta-sample of 2300 AGN with ~ 50% redshift com- 
pleteness and find a shallow increase of the type-2 AGN fraction 
proportional to (1 + z) 0A . The differences between the various 
results are very likely due to the different sample sizes and the 
way, in which the redshift incompleteness, the different flux lim- 
its and the AGN classification have been taken into account in 
the analysis. 

The analysis in this paper shows a significant trend of an 
increasing absorbed fraction with redshift. Interestingly, at first 



sight, the observed fraction of absorbed sources, after correction 
for the identification incompleteness has a redshift dependence 
very similar to the prediction of the no evolution GCH07 model 
(see Figure [7] right). This, however, turns out to be a "conspir- 
acy" due to a combination of the shape of the luminosity- and 
redshift-dependence of the absorbed fraction. As can be seen 
in Figure [3] at higher redshifts the median of the sample shifts 
to higher luminosities. In the case of the almost linear decrease 
of the absorbed fraction with luminosity observed in Figure [7j 
(left) one would expect a significant decrease of the observed 
absorbed fraction if there was no redshift evolution (see also the 
model predictions in Figure 2 of the Treister & Urry 2006 paper). 
This has to be compensated by an increase of the absorbed frac- 
tion with redshift if an almost constant behaviour like the one in 
Figure [7j (right) is observed. Alternatively, if at high luminosi- 
ties the absorbed fraction is constant, as assumed in the GCH07 
model, an almost constant trend with redshift is expected, as- 
suming no evolution. 

However, disentangling the data into separate redshift and lu- 
minosity intervals resolves the degeneracy of the analyses based 
on one dimensional analyses. Figure|8]confirms with high statis- 
tical accuracy that the decreasing trend of the absorbed fraction 
with luminosity is present in all redshift intervals. A similar re- 
sult has already been shown in broader redshift and luminosity 
bins by Barger et al. (2005), however, their results are affected by 
the systematic differences discussed in subsection l7.ll However, 
the Figure[8]also shows an increase of the absorbed fraction with 
redshift. If the data in each redshift slice are fit by a simple linear 
relation, the normalization of this relation, here assumed for a lu- 
minosity log(Lx) = 43.75 in the middle of the observed range, 
clearly shows an increase with redshift up to z~2, where this 
trend seems to saturate. Figure [9] shows that this trend is highly 
significant and can be fit by a power law evolution proporional 
to (l+z) ff , with a = 0.62 + 0.11. This evolution is thus substan- 
tially faster than the exponent a = 0.4 ±0.1 originally found 
by Treister & Urry (2006). Several reasons can be envisaged for 
this difference: First, these authors have fit the power law trend 
over the full redshift range, because the flattening of the curve 
was not apparent in their data. In addition they used the BLAGN 
classification only, which thends to over-estimate the absorbed 
fraction at low redshifts (see Figure [TOb. Finally, the inclusion 
of crude redshifts in the analysis presented here tends to recover 
more absorbed objects at higher redshifts and higher luminosi- 
ties. 

Interestingly, the increase of the normalization of the ab- 
sorbed fraction relation with redshift is very similar to the ob- 
served fractions given in figure 6b of the the La Franca et al. 
(2005) paper based mainly on the HELLAS2XMM sample (see 
also Figure |5). This could, however, be somewhat of a chance 
coincidence, because their data shows the observed ratios in- 
tegrated over a broad luminosity range (more like Figure |7j 
(right) in this paper), and has its own systematic selection ef- 
fects. Nevertheless, the trend originally claimed by La Franca et 
al. d2005 ) is fully confirmed here. 

It is also interesting to note that the saturated absorption evo- 
lution found in this study shows a redshift dependence very sim- 
ilar to the space density evolution of the overall AGN popula- 
tion with an increase from zero redshift to a peak redshift, where 
the evolution is saturating. However, the absorption evolution is 
much shallower than the AGN density evolution. 
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7.4. Consequences for the broader context of AGN and 
galaxy co-evolution 

The cosmological evolution of AGN in the X-ray, optical and 
radio wavebends can be described by a luminosity-dependent 
density evolution model, in which the peak of the AGN space 
density shifts to lower luminosities towards lower redshifts (see 
above and Hopkins et al. 120071 ). This evolutionary behaviour is 
very similar to the "cosmic downsizing" observed in the bolo- 
metric luminosity of the normal galaxy population (Cowie et al. 
1996 ) and indicates that star formation and nuclear activity in 
galaxies go hand in hand. This kind of "anti-hierarchical" Black 
Hole growth scenario is not predicted by most semi-analytic 
AGN evolution models based on galaxy merger scenarios (e.g. 
K auffma nn & Haehnelt 2000, Wyithe & Loeb [200l Croton et 
al. 120031 Menci et al. 120081 Rhook & Haehnelt[2008]). Typically 
these models are able to explain the space density evolution of 
high-luminosity AGN, peaking at redshifts z>2, but have dif- 
ficulties to match the behaviour at lower luminosities. Some 
authors (e.g. Treister et al. 120041 Menci et al. 12008) attribute 
these differences to systematic selection biases against the iden- 
tification of distant obscured, low-luminosity AGN, in partic- 
ular if there is a strong evolution of absorption with redshift. 
The results of the analysis presented in this paper, although in- 
deed confirming a significant evolution of obscuration, argue 
against this conjecture, because the evolution observed across 
the whole redshift range is too shallow (at most a factor of 2) 
to explain the dramatic luminosity-dependent density evolution 
effect. Therefore other ingredients are required to explain the ob- 
served AGN downsizing. 
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Fig. 11. Schematic arrows showing galaxies moving in the 
colour-mass diagram under different evolutionary assumptions 
(following Faber et al.. 120071 ). It is conventionally assumed that 
galaxies on the "red sequence" are formed when two "blue 
cloud" galaxies merge with each other and the feedback from 
the growing central black hole is quenching star formation, leav- 
ing a "dead" and red bulge-dominated galaxy with a dormant su- 
permassive black hole. Host galaxies of intermediate-luminosity 
active galactic nuclei are often found in the region between the 
blue cloud and red sequence, the so called "green valley". The 
scenario described here tries to explain this by a "re-juvenation" 
of bulge-dominated galaxies through the accretion of fresh gas 
from their environment. 



Important inputs to the solution of this puzzle may come 
from the study of colours and morphology of the host galaxies 
of moderate-luminosity AGN at low to intermediate redshifts. 
The SDSS in combination with data from the Galaxy Evolution 
Explorer (GALEX) satellite are of help here. Kauffmann et al. 
(2007) analyzed a volume-limited sample of massive bulge- 
dominated galaxies with data from both the SDSS and GALEX 
at redshifts 0.03<z<0.07. The GALEX NUV data are a very 
sensitive indicator for low levels of star formation in these sys- 
tems. These authors find that practically all AGN in their sam- 
ple have rather blue NUV-r colours and the UV excess light 
is almost always associated with an extended blue disk, while 
galaxies with red outer regions almost never have a young bulge 
or a strong AGN. They suggest a scenario in which the gas of 
the outer disk is accreted onto the bulge-dominated galaxy over 
cosmic time and provides the mass reservoir to trigger sporadic 
AGN accretion and star formation through gravitational instabil- 
ities (see also Cavaliere & Vittorini 2000; Hammer et al. 2005 ). 
Interestingly, Schiminovich et al. (2007 ) show that in the NUV-r 
versus mass diagram, the local SDSS AGN are preferentially at 
high masses, but at colours in between the blue cloud and the red 
sequence, the so called "green valley" (see also Figure ITTb. 

Similar results are found at higher redshifts from optical 
studies in deep and wide Chandra and XMM-Newton survey 
fields. Nandra et al. (120071 ) discuss the colour-luminosity dia- 
gram for AGN selected from the Chandra survey of the Extended 
Groth Strip in the redshift range 0.6<z<1.4. They also find the 
AGN host galaxies in a distinct region of the colour-magnitude 
diagram, often in the "green valley". They interprete this in the 
context of star formation quenching in massive galaxies, result- 
ing in a migration from the blue cloud to the red sequence. 
Silverman et al. (2008) present an analysis of 109 moderate- 
luminosity (log Lx <42) AGN in the Extended Chandra Deep 
Field-South survey (Lehmer at al. 2000) with redshifts z=0.4- 
1.1 and find a strong dependence of AGN host colours on the 
environment. Compared to all galaxies in the field they find that 
the fraction of galaxies hosting an AGN peaks in the "green val- 
ley", in particular due to enhanced AGN activity in two narrow 
redshift spikes around z=0.63 and z=0.73 (see Gilli et al. 2005 ). 
They also find that AGN host galaxies in this redshift range typ- 
ically have hybrid morpholgies between pure bulges and disks. 
They try to explain their findings in an AGN merger scenario, 
but remark that the merger timescales are typically too short to 
explain the abundance of green valley galaxies. 

Putting these findings together, one arrives at a scenario pos- 
sibly explaining the AGN cosmic downsizing as well as their 
colour and morphological evolution (see Figure [TTll. In order to 
produce an active galactic nucleus, we need two ingredients: a 
central supermassive black hole and efficient gas accretion. The 
fuel can either be provided by a major merger driving gas into 
the center and simultaneously leading to an Eddington-limited 
growth of the black hole(s) (see e.g. Di Matteo et al. 120051 Li 
et al. 120061 ). In this case a bulge galaxy is formed and the feed- 
back of the growing black hole can quench star formation. The 
result is a red bulge-dominated galaxy with a dormant supermas- 
sive black hole. The merger time scale is, however, rather short 
(~ 10 8 years), so that only a very small fraction of the galaxy 
population can be found in a transient state between the blue 
cloud and the red sequence. It is likely that the space density evo- 
lution of the most luminous QSOs, peaking in the redshift range 
z=2-3, is a direct consequence of the merger history of massive 
halos. The results above show that there is, however, also another 
source of gas supply, namely the accretion of fresh gas from the 
surrounding of a galaxy. In this way, a bulge-dominated galaxy 
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with a central black hole which was formed during a previous 
merger event can be "re-juvenated", leading to a disk surround- 
ing the bulge. In the colour-mass diagram such a galaxy is there- 
fore coming down from the red sequence into the green valley 
(see Figure ITTb and the galaxy assumes the hybrid bulge/disk 
morphology observed for the SDSS AGN and many of the X- 
ray selected intermediate-redshift AGN. The black hole can be 
fed through episodic star formation and accretion events, likely 
produced by gravitational instabilities in the outer disk. Since 
the black hole and the bulge have already been formed previ- 
ously and the host galaxy is already relatively massive, the ac- 
cretion typically happens at a rather low Eddington ratio, sub- 
stantially smaller than in the merger case (see e.g. Marconi et al. 
120041 Merloni & Heinz l2008l l. This mode of black hole growth is 
therefore slow enough that it can be observed in a large fraction 
of all galaxies, providing a natural explanation for the late evo- 
lution in the ubiquitous lower-luminosity AGN. This scenario is 
still rather descriptive and crude, but hopefully can be quanti- 
fied by including the bi-modal gas supply into semi-analytical 
evolution models in the future. 



8. Conclusions 

1 . A meta-sample of hard X-ray selected AGN has been com- 
piled from the literature, providing an unprecedented com- 
bination of statistical quality and spectroscopic as well as 
photometric redshift completeness. Systematic differences in 
the selection of this rather heterogeneous set of samples were 
taken care of as far as possible and should not significantly 
bias the results. 

2. The small, but significant redshift incompleteness could be 
corrected for using the correlation between luminosity and 
X-ray to optical flux ratio observed in non-broad-line AGN. 
A comparison between the analysis with and without the in- 
clusion of crude redshifts carried throughout the paper shows 
that the crude redshifts tend to enhance the high-redshift 
high-luminosity bins, but the main results do not depend 
strongly on this choice. 

3. The X-ray selected AGN could be classified using both opti- 
cal spectroscopy and X-ray hardness ratios. This classifica- 
tion scheme appears more robust than either optical or X-ray 
classification alone. 

4. A strong decrease of the fraction of absorbed AGN was 
found with redshift. This trend confirms similar results ob- 
tained previously in the X-ray, but also optical and MIR 
bands. As expected, systematic selection effects are present 
between different bands, but do not dominate the results. 
This signals a break-down of the strong unified AGN model 
and indicates that high-luminosity AGN can clean out their 
environment. 

5. The same decreasing trend with luminosity is found in all 
individual redshift shells up to z~3. However, the data in- 
dicate a significant evolution of the absorbed AGN fraction, 
which increases by about a factor of 2 up to z~2 and remains 
approximately constant thereafter. 

6. This rather shallow evolution implies that the luminosity- 
dependent AGN density evolution observed in the X-ray, 
optical and radio bands is a real property of the parent popu- 
lation and not caused by systematic selection effects. 

7. A new scenario is proposed, in which the cosmic down- 
sizing of the AGN population is due to two different fuelling 
mechanisms, on one hand the efficient growth of luminous 
QSOs by galaxy mergers early in the universe, on the other 



hand the re-juvenation of pre-formed bulges and black holes 
by slow gas accretion over cosmic time. 
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